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ABSTRACT 

Context. The Spitzer Space Telescope has detected a powerful (Lh, ~ 10 41 erg s" 1 ) mid-infrared H2 emission towards the galaxy-wide 
collision in the Stephan's Quintet (SQ) galaxy group. This discovery was followed by the detection of more distant H2-luminous 
extragalactic sources, with almost no spectroscopic signatures of star formation. These observations place molecular gas in a new 
context where one has to describe its role as a cooling agent of energetic phases of galaxy evolution. 

Aims. The SQ postshock medium is observed to be multiphase, with H 2 gas coexisting with a hot (~ 5x 10 6 K), X-ray emitting plasma. 
The surface brightness of H2 lines exceeds that of the X-rays and the 0-0 S( 1 ) IT linewidth is ~ 900 km s~' , of the order of the collision 
velocity. These observations raise three questions we propose to answer: ( i) Why H 2 is present in the postshock gas ? ( ii) How can we 
account for the IT excitation ? ( Hi) Why is IT a dominant coolant ? 

Methods. We consider the collision of two flows of multiphase dusty gas. Our model quantifies the gas cooling, dust destruction, H2 
formation and excitation in the postshock medium. 

Results, (i) The shock velocity, the post-shock temperature and the gas cooling timescale depend on the preshock gas density. The 
collision velocity is the shock velocity in the low density volume filling intercloud gas. This produces a ~ 5 x 10 6 K, dust-free, 
X-ray emitting plasma. The shock velocity is smaller in clouds. We show that gas heated to temperatures less than 10 6 K cools, keeps 
its dust content and becomes IT within the SQ collision age (~ 5 x 10 6 years), (ii) Since the bulk kinetic energy of the H2 gas is the 
dominant energy reservoir, we consider that the H 2 emission is powered by the dissipation of kinetic turbulent energy. We model this 
dissipation with non-dissociative MHD shocks and show that the H2 excitation can be reproduced by a combination of low velocities 
(in the range 5-20 km s _1 ) shocks within dense (n H > 10 3 cm~ 3 ) H 2 gas. (Hi) An efficient transfer of the bulk kinetic energy to 
turbulent motions of much lower velocities within molecular gas is required to make H 2 a dominant coolant of the postshock gas. We 
argue that this transfer is mediated by the dynamical interaction between gas phases and the thermal instability of the cooling gas. 
We quantify the mass and energy cycling between gas phases required to balance the dissipation of energy through the H 2 emission 
lines. 

Conclusions. This study provides a physical framework to interpret H2 emission from H2-luminous galaxies. It highlights the role that 
H2 formation and cooling play in dissipating mechanical energy released in galaxy collisions. This physical framework is of general 
relevance for the interpretation of observational signatures, in particular H 2 emission, of mechanical energy dissipation in multiphase 
gas. 

Key words. ISM: general - ISM: dust, extinction - ISM: molecules - Atomic processes - Molecular processes - Shock waves - 
Galaxies: evolution - Galaxies: interactions 



1. Introduction 

Spitzer Infra-Red Spectrograph (IRS) observations of the 
Stephan's Quintet (hereafter SQ) galaxy-wide collision led to 
the unexpected detection of extremely bright mid-IR H2 ro- 
tational line emission from warm (~ 10 2 - 10 3 K) molecu- 
lar gas (Appleton et al., 2006). This is the first time an almost 
pure H2 line spectrum has been seen in an extragalactic ob- 
ject. This result is surprising since H2 is coexisting with a hot 
X-ray emitting plasma and almost no spectroscopic signature 
(dust or ionized gas lines) of star formation is associated with 
the H2 emission. This is unlike what is observed in star form- 
ing galaxies where the H2 lines are much weaker than the mid- 
IR dust features (Rigopoulou et al., 2002; Higdon et al., 2006; 
Roussel et al., 2007). 
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The detection of H2 emission towards the SQ shock was 
quickly followed by the discovery of a class of ^-luminous 
galaxies which show high equivalent width H2 mid-IR lines. The 
first H^-luminous galaxy, NGC 6240, was identified from ground 
based near-IR H 2 spectroscopy by Joseph et al. (1984), although 
this galaxy also exhibits a very strong IR continuum. Many more 
of H2 bright galaxies are being found with Spitzer (Egami et al., 
2006; Ogle et al., 2007), often with very weak continuua sug- 
gesting that star formation is not the source of H2 excitation. 
These Kh-luminous galaxies open a new perspective on the role 
of H2 as a cold gas coolant, on the relation between molecular 
gas and star formation, and on the energetics of galaxy forma- 
tion, which has yet to be explored. 

SQ is a nearby H2-luminous source where observations pro- 
vide an unambigious link between the origin of the H2 emission, 
and a large-scale high-speed shock. The SQ galaxy-wide shock 
is created by an intruding galaxy colliding with a tidal tail at a 
relative velocity of ~ 1 000 km s _1 . Evidence for a group-wide 
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shock comes from observations of X-rays from the hot postshock 
gas in the ridge (Trinchieri et al., 2003, 2005; O'Sullivan et al., 
2008), strong radio synchrotron emission from the radio emit- 
ting plasma (Sulentic et al., 2001) and shocked-gas excitation 
diagnostics from optical emission lines (Xu et al., 2003). Spitzer 
observations show that this gas also contains molecular hydro- 
gen and that it is turbulent with an H2 linewidth of 870 km s _1 . 
The H2 surface brightness is larger than the X-ray emission from 
the same region, thus the H2 line emission dominates over X-ray 
cooling in the shock. As such, it plays a major role in the energy 
dissipation and evolution of the postshock gas. 

These observations raise three questions we propose to an- 
swer: ( i) Why is there H2 in the postshock gas ? ( ii) How can 
we account for the H2 excitation ? (Hi) Why is Ho a dominant 
coolant ? We introduce these three questions, which structure 
this paper. 

( i) Hi formation. The detection of large quantities of warm 
molecular gas in the SQ shock, coexisting with a X-ray emitting 
plasma, is a surprising result. Appleton et al. (2006) invoked an 
oblique shock geometry, which would reduce the shock strength. 
However, with this hypothesis, it is hard to explain why the tem- 
perature of the postshock plasma is so high (5 x 10 6 K, which 
contrains the shock velocity to be ~ 600 km s ). Therefore the 
presence of IT is likely to be related to the multiphase, cloudy 
structure of the preshock gas. 

One possibility is that molecular clouds were present in the 
preshock gas. Even if the transmitted shock into the cloud is dis- 
sociative, H2 molecules may reform in the postshock medium 
(Hollenbach & McKee, 1979). An alternative possibility is that 
H2 forms out of preshock Hi clouds. Appleton et al. (2006) 
proposed that a large-scale shock overruns a clumpy preshock 
medium and that the H2 would form in the clouds that experi- 
ence slower shocks. In this paper we quantify this last scenario 
by considering the collision between two inhomogeneous gas 
flows, one being associated with the tidal tail and the other as- 
sociated with the interstellar medium (hereafter ISM) of the in- 
truding galaxy. 

(ii) Ho excitation. Several excitation mecanisms may ac- 
count for high-equivalent width H2 line emission. The excita- 
tion by X-ray photons was quantified by a number of authors 
(e.g. Tine et al., 1997; Dalgarno et al., 1999). The energy con- 
version from X-ray flux to H2 emission is at most 10% for a 
cloud that is optically thick to X-ray photons. The absorbed frac- 
tion of X-ray photons may be even smaller if the postshock H2 
surface filling factor is smaller than 1 . Chandra and XMM obser- 
vations (0.2 - 3 keV) show that, within the region where Spitzer 
detected H2 emission, the H2 to X-ray luminosity ratio is ~ 5 
(O'Sullivan et al., 2008; Trinchieri et al., 2005). Therefore, the 
excitation of H2 by 0.2 - 3 keV X-ray photons cannot be the 
dominant process. 

H2 excitation may also be produced by cosmic ray ionization 
(Ferland et al., 2008). However radio continuum observations of 
SQ show that the combined cosmic ray plus magnetic energy is 
not the dominant energy reservoir. 

SQ observations suggest that only a fraction of the colli- 
sion energy is used to heat the hot plasma. Most of this en- 
ergy is observed to be kinetic energy of the H2 gas. Therefore 
we consider that the H2 emission is most likely to be pow- 
ered by the dissipation of the kinetic energy of the gas. This 
excitation mechanism has been extensively discussed for the 
Galactic interstellar medium. It has been proposed to account 
for H2 emission from solar neighbourhood clouds (Gry et al., 
2002; Nehme et al., 2008), from the diffuse interstellar medium 



in the inner Galaxy (Falgarone et al., 2005) and from clouds in 
the Galactic center (Rodriguez-Fernandez et al., 2001). 

( Hi) Hi cooling. For H2 to be a dominant coolant of the SQ 
postshock gas, the kinetic energy of the collision in the center of 
mass rest-frame has to be efficiently transfered to the molecular 
gas. In the collision, the low density volume filling gas is decel- 
erated but the clouds keep their preshock momentum. The clouds 
which move at high velocity with respect to the background 
plasma are subject to drag. This drag corresponds to an exchange 
of momentum and energy with the background plasma. This 
transfer of energy between clouds and hot plasma has been dis- 
cussed in diverse astrophysical contexts, in particular infalling 
clouds in galaxy clusters cooling flows (Pope et al., 2008) and 
galactic halos (Murray & Lin, 2004), as well as cometary clouds 
in the local interstellar medium (Nehme et al., 2008) and the 
Helix planetary nebula tails (Dyson et al., 2006). 

The dynamical evolution of the multiphase interstellar 
medium (ISM) in galaxies has been extensively investigated with 
numerical simulations, within the context of the injection of 
mechanical energy by star formation, in particular supernovae 
explosions (e.g. de Avillez & Breitschwerdt, 2005a; Dib et al., 
2006). The injection of energy by supernovae is shown to be 
able to maintain a fragmented, multiphase and turbulent ISM. 
Numerical simulations on smaller scales show that the dy- 
namical interaction between gas phases and the thermal in- 
stability feed turbulence within clouds (Sutherland et al., 2003; 
Audit & Hennebelle, 2005; Nakamura et al., 2006). The SQ col- 
lision releases ~ 10 56 erg of mechanical energy in a collision 
timescale of a few million years. This is commensurate with 
the amount of mechanical energy injected by supernovae in the 
Milky Way (MW) over the same timescale. The main difference 
is in the mass of gas, which is two orders of magnitude smaller 
in the SQ shock than in the MW molecular ring. None of the 
previous simulations apply to the context of the SQ collision, 
but they provide relevant physics to our problem. In this paper, 
based on this past work, we discuss how the dynamical inter- 
action between the molecular gas and the background plasma 
may sustain turbulence within molecular clouds at the required 
level to balance the dissipation of energy through the H2 emis- 
sion lines. 

The structure of this paper is based on the previous 3 ques- 
tions and is organised as follows: in § 2, we gather the observa- 
tional data that set the luminosity, mass and energy budget of the 
SQ collision. In § 3, the gas cooling, dust destruction, and H2 
formation in the postshock gas are quantified. In § 4 the dissi- 
pation of kinetic energy within the molecular gas is modeled, in 
order to account for the H2 emission in SQ. Section 5 presents a 
view at the dynamical interaction between ISM phases that arises 
from our interpretation of the data. In § 6 we discuss why H2 is 
such a dominant coolant in the SQ ridge. Section 7 discusses 
open questions including future observations. Our conclusions 
are presented in § 8. 

2. Observations of Stephan's Quintet 

This section introduces the observational data that set the lu- 
minosity, mass and energy budgets of the SQ shock. Relevant 
observational numbers and corresponding references are gath- 
ered in Table 1. Note that all the quantities are scaled to the 
aperture J{ used by Appleton et al. (2006) to measure H2 lu- 
minosities : J{ = 11.3 x 4.7 arcsec 2 , which corresponds to 
5.2x2.1 kpc 2 . The equivalent volume associated with this area is 
'Vw = 5.2 x 2.1 x l z kpc 3 , where I- [kpc] is the dimension along 
the line-of-sight. For the X-ray emitting gas, O'Sullivan et al. 
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(2008) assume a cylindrical geometry for the shock, which gives 
L ~ 4.5 kpc. This dimension may be smaller for the H2 emit- 
ting gas, which could be concentrated where the gas density is 
largest. To take this into account, we use a reference value of 
L = 2 kpc. Since this number is uncertain, the explicit depen- 
dence of the physical quantities on this dimension are written 
out in the equations. 



2. 1 . Astrophysical context 

Since its discovery (Stephan, 1877), the Stephan's Quintet 
(94 Mpc) multigalaxy compact group has been extensively stud- 
ied, although one member of the original quintet (NGC 7320) 
is now presumed to be a foreground dwarf and will not be dis- 
cussed further here. Later studies have restored the "quintet" sta- 
tus by the discovery of a more distant member. The inner group 
galaxies are shown in the left-hand side of Fig. 1. 

Hi observations exhibit a large stream created by tidal in- 
teractions between NGC 7320c, and NGC 7319 (Moles et al., 
1997; Sulentic et al., 2001). It is postulated that another galaxy, 
NGC 7318b (hereafter called the "intruder"), is falling into the 
group with a very high relative velocity shocking a 25 kpc seg- 
ment of that tidal tail (Xu et al., 2003). The observed differ- 
ence in radial velocities between the tidal tail and the intruder 
is ~ 1 000 km s _1 . If this picture of two colliding gas flows is 
right, we expect two shocks: one propagating forward into the 
SQ intergalactic medium, and an other reverse shock driven into 
the intruding galaxy. 

Over the shock region, no Hi gas is detected, but optical 
line emission from Hn gas is observed at the Hi tidal tail veloc- 
ity. X-ray (Trinchieri et al., 2003, 2005; O'Sullivan et al., 2008) 
observations resolve shock-heated plasma spatially associated 
with the optical line and 21 cm radio continuum (Sulentic et al., 
2001). The width of this emission is ~ 5 kpc, which is commen- 
surate with the width of the H a emission. Dividing this shock 
width by the collision velocity (assumed to be 1 000 km s _1 ), we 
get a collision age of t co \\ ~ 5 x 10 6 yr. 

Spectral energy distribution fitting of Chandra and XMM 
observations of SQ show that the temperature of the hot 
plasma is T x ~ 5 x 10 6 K (Trinchieri et al., 2003, 2005; 
O'Sullivan et al., 2008). The postshock temperature is given by 
(see e.g. Draine & McKee, 1993): 



_ 2(7-1) n 
ps (y + D 2 k B 



^V 2 -5xl0 6 ( ^ r ) K, 

k B \600kms- 1 ; 



(1) 



where V s is the velocity of the shock wave, fi the mean par- 
ticle mass (equals to 10~ 24 g for a fully ionized gas), k B the 
Boltzmann constant and y = 5/3. Therefore, the postshock tem- 
perature of the hot plasma allows to estimate a shock velocity 
of ~ 600 km s . This shock velocity is consistent with the gas 
velocity in the center of mass frame, which would be approxi- 
mately half of the observed relative velocity if the mass of the 
gas reservoirs on the intruder and SQ tail sides are commensu- 
rate. 



2.2. Mass budgets 

Observations show that the preshock and postshock gas are mul- 
tiphase. We combine H2, Hi, Hn and X-ray gas luminosities to 
estimate gas masses. 



2.2.1. Preshock gas 

The Hi preshock gas is contained on both the SQ tidal tail side 
(at a velocity of ~ 6 700 km s _1 ) and the intruder galaxy side 
(~ 5 700 km s _1 ). On both sides, this gas is seen outside the 
shock area (Sulentic et al., 2001; Williams et al., 2002). On the 
SQ intra-group side, the tidal tail Hi column densities are in the 
range Nn = 1 - 3 X 10 20 cm 2 . These two values bracket the 
range of possible preshock column densities. Neutral hydrogen 
observations of SQ show Hi at the intruder velocity to the South 
West (SW) of the shock area. The Hi column density of the SW 
feature and that of the tidal tail are comparable. By multiplying 
the column densities by the area Si, we derive a total mass of Hi 
preshock gas in the aperture Si of 2 - 6 x 10 7 M . 

Deep Chandra and XMM-Newton observations show a dif- 
fuse "halo" of soft X-ray emission from an extended re- 
gion around the central shock region (Trinchieri et al., 2005; 
O'Sullivan et al., 2008). We consider that this emission is trac- 
ing the preshock plasma, which may have been created by 
a previous collision with NGC 7320c (Sulentic et al., 2001; 
Trinchieri et al., 2003). The mass of hot gas, within the equiv- 
alent volume Vji associated with the aperture Si, is derived from 
the X-ray luminosity and the temperature Tx- 
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1/2 (3 x 10" 



23 \!/ 2 



A(7x) 



I, 



2 kpc 
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where Lx [erg s~'] is the X-ray luminosity, A(7x) [erg s _1 cm 3 ] 
the gas cooling efficiency at the plasma temperature Tx- The X- 
ray luminosity in the halo is derived from Trinchieri et al. (2005) 
and is 2.5 x 10 39 erg s -1 . We assume a solar metallicity (Xu et al., 
2003) and a cooling efficiency of A(7x) = 3 x 10~ 23 erg s _1 cm 3 
at Tx = 5 x 10 6 K (see Appendix B.l and Gnat & Sternberg, 
2007, for more details). 

2.2.2. Postshock gas 

The galaxy-wide collision creates a multi-phase interstellar 
medium where molecular and Hn gas are embedded in a hot X- 
ray emitting plasma. 

No postshock Hi gas is detected in the shocked region where 
H2 was observed. We derive an upper limit of 5 x 10 6 M based 
on the lowest Hi contour close to the shock. A two-temperature 
fit of the H2 excitation diagram gives T = 185 K and T = 675 K 
for the two gas components. Most of the gas mass is contained 
in the lower temperature component, which represents a mass of 
warm H 2 of 3 x 10 7 M o (Apple ton et al., 2006) within the aper- 
ture Si. Published CO observations do not constrain the total H2 
mass in the preshock and postshock gas. Given the uncertainties 
on the preshock gas masses, one cannot infer a precise mass of 
cold molecular gas. 

The mass of hot plasma is derived from the O'Sullivan et al. 
(2008) Chandra observations in the region observed by Spitzer. 
The absorption-corrected X-ray luminosity is 1.2 x 10 40 erg s _1 
within the aperture Si. Using equation 2, this corresponds to a 
mass of 1.4x 10 7 M . This is 2-3 times the mass of the preshock 
hot plasma. O'Sullivan et al. (2008) report variations of the av- 
erage proton density in the range n H = 1-2 - 2.3 x 10 -2 cnT 3 . 
The density peaks close to the area Si. We adopt a hot plasma 
density of «h - 0.02 cirT 3 where H2 has been detected. This is 
in agreement with the value obtained by Trinchieri et al. (2005). 
The postshock gas pressure is given by 



^=2.2xl0 5 °- 6XOTH 
k B fi 



(—?»—){ 
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5 x 10 6 K 
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Preshock Gas 


Postshock Gas 


Density n H [cm -3 ] 
Temperature T [K] 
Pressure P/ku [ K cm" 3 ] 
N H [cm- 2 ] 


Hot Plasma" Hi c H 2 
3 x 10 20 


Hot Plasma" Hn b Hi c WarmH 2 d Cold H 2 
0.02 

5xl0 6 185 < 10 2 
2x 10 5 

1.4 xlO 19 <5.8xl0 19 2xl0 20 


Mass Budget [M G ] 


Hot Plasma" Hi c H 2 
6 x 10 6 2 - 6 x 10 7 


Hot Plasma" Hn b Hi c Warm H 2 d Cold H 2 
1.4 xlO 7 1.2 xlO 6 <5xl0 6 3 xl0 7 


Energy Budget [erg] 


Thermal Bulk Kinetic 
(Plasma in halo) (Shocked H 1 gas) 
1.7 x 10 55 1 - 3 x 10 56 


Thermal Bulk Kinetic e 
(X-ray emitting plasma) (Turbulent H 2 ) 
4xl0 55 > 1.2 xlO 56 


Flux [W nr 2 ] 


X-rays 
2.3 x 1(T 18 


X-rays Ka Oi H 2 d 
1.1 xlO -17 5.1xlO~ 18 3.6 x 10~ 18 5.5xl0~ 17 


Luminosity f [erg s ] 


2.5 x 10 39 


1.2 xlO 40 5.4 xlO 39 3.8 xlO 39 5.8 x 10 40 



Table 1. Mass, energy and luminosity budgets in front and behind the Stephan's Quintet galaxy-wide shock, directly inferred 
from observations. All the numbers are scaled to our aperture = 5.2 x 2.1 kpc 2 ) where Spitzer observations were performed. 
The preshock gas is mainly the Hi gas contained in the tidal tail of NGC 7319. After the shock, the mass is distributed between 
the hot X-ray emitting plasma and the H2 gas (see text for details). Before the shock, most of the energy available is the kinetic 
energy of the Hi gas that will be shocked. The shock splits the energy budget in two parts: thermal and kinetic energy. The former 
is stored in the hot plasma whereas the latter goes in turbulent motions that heat the H2 gas. Observations show that mechanical 
energy is the dominant energy reservoir. 



" XMM-Newton observations of the extended X-ray emission in the shock and the tidal tail (Trinchieri et al., 2005). 
b Ha and Oi optical line observations by Xu et al. (2003). 

Based on extrapolation of Hi observations in the tidal tail by Sulentic et al. (2001); Williams et al. (2002). 

d From Spitzer IRS observations. The H 2 line flux is summed over the S(0) to S(5) lines (Appleton et al., 2006). 

e This bulk kinetic energy is a lower limit because an unknown mass of turbulent cold molecular gas can contribute to this energy. 

t Luminosities are indicated assuming a distance to SQ of 94 Mpc. 



In the tenuous gas this pressure remains almost constant during 
the collision time scale since the plasma does not have the time 
to cool nor to expand significantly. The isobaric gas cooling time 
scale can be estimated as the ratio of the internal energy of the 
gas to the cooling rate: 

5/2fc B r x _ in7 / T x \ 1M /0.02cm- 3 \ 

fcooi= ^7-^9x10' - — — — yr.(4) 

« H A(r x ) \5xl0 6 K/ \ « H / 

This cooling time is significantly longer than the age of the 
shock. In the right hand side of Eq. 4, we use the power-law fit of 
the Z = 1 cooling efficiency given in Gnat & Sternberg (2007) 
(consistent with Sutherland & Dopita, 1993). If the metallicity 
of the pre-shocked tidal filament was sub-solar, this would in- 
crease the cooling time of the gas still further. O' Sullivan et al. 
(2008) found a best fit to their X-ray spectrum for Z = 0.3 Z . 

The mass of the ionized Hn gas is derived from the emission 
measure of H a observations: 

Mh^1.2xhW J(^n— K (5) 

\5.1 x 10- 18 /\10 4 K/ \ « H / 

where Tna [W m 2 ] is the flux of the H a line scaled to the aper- 
ture J{ and « H the gas proton density. For SQ, the observed H a 
flux is Thv = 5.1 x 10~ 18 W nr 2 (Xuet al., 2003). At the SQ 
postshock pressure, the gas density at 10 4 K is «h - 10 cirT 3 . 

2.3. Extinction and UV field 

Xu et al. (2003) have interpreted their non-detection of the Hg 
line as an indication of significant extinction (E(Hp - H ff ) > 1 .2, 



i.e. A(H a ) > 2.4). However, more sensitive spectroscopy 1 by R- 
A. Due et al. (private communication) show that the Hg line is 
detected over the region observed with IRS at two velocity com- 
ponents (~ 5 700 and 6 300 km s _1 , corresponding to the intruder 
and the intra-group gas velocities). These observations show that 
the extinction is much higher (Ay = 2.5) for the low-velocity gas 
component than for the high-velocity one. There are also large 
spatial variations of the extinction value for the intruder gas ve- 
locity close to this position, suggesting dumpiness of the gas 
within the disturbed spiral arm of NGC 7318b. An extinction 
of Ay = 0.3 - 0.9 is derived from the Balmer decrement for 
the intra-group tidal tail velocity, with much lower spatial vari- 
ations. After subtraction of the foreground galactic extinction ( 
Ay = 0.24), this gives a SQ ridge extinction Ay in the range 
0.1-0.7. 

A low extinction for the H2-rich, intra-group gas is also 
consistent with measured UV and far-IR dust surface bright- 
ness observations. Using the published value of FUV extinction 
(Xuetal., 2005) and the Weingartner & Draine (2001) extinc- 
tion curve, we derive Ay = 0.29, which is in agreement with the 
new B aimer-decrement study. GALEX observations show that 
the radiation field is G ~ 1 - 2 Go in the shocked region, where 
Go is the interstellar radiation field (ISRF) in Habing units 2 . 
We obtain the higher value using the flux listed in Table 1 of 
Xu et al. (2005). This flux is measured over an area much larger 
than our aperture It provides an upper limit on the UV field 
in the shock because this area includes star forming regions in 



1 observations made at the Cerro Tololo Inter- American Observatory 
(CTIO) 

2 Go = 2.3 x 10~ 3 erg s" 1 enr 2 at 1 530 A, Habing (1968) 
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the spiral arm of NGC 7318b. ISO observations show a dif- 
fuse 100 /Yin emission at 2.2 MJy sr~' level in the shocked re- 
gion (Xu et al., 2003). Using the Galactic model of Draine & Li 
(2007) for G = 2 Go, the 100 /mi FIR brightness corresponds to 
a gas column density of 1.5 x 10 20 cirT 2 . For G = Go, we obtain 
N H = 3.3 x 10 20 cirT 2 , which corresponds to A v = 0.18 for the 
Solar neighbourhood value of the Av/Nh ratio. Since this value 
is low, no extinction correction is applied. 

2.4. Energy reservoirs 

The bulk kinetic energy of the pre- and postshock gas are 
the dominant terms of the energy budget of the SQ collision 
(Table 1). Assuming that the preshock Hi mass is moving at 
V s = 600 km s _1 in the center of mass rest frame, the preshock 
gas bulk kinetic energy is 1 - 3 x 10 56 erg. Note that this is a 
lower limit since a possible contribution from preshock molecu- 
lar gas is not included. The thermal energy of the X-ray plasma 
represents ~ 10 % of the preshock kinetic energy. 

The shock distributes the preshock bulk kinetic energy into 
postshock thermal energy and bulk kinetic energy. The Spitzer 
IRS observations of SQ (Appleton et al., 2006) show that the 0- 
S(l) 17 /mi H 2 line is resolved with a width of 870 + 60 km s -1 
(FWHM), comparable to the collision velocity. We assume that 
the broad resolved line represents H2 emission with a wide range 
of velocities, by analogy with the Oi [6 300 A] emission seen in 
the same region by Xu et al. (2003) (see Appleton et al. (2006) 
for discussion). Then the results imply a substantial kinetic en- 
ergy carried by the postshock gas. The kinetic energy of the H2 
gas is estimated by combining the line width with the warm H2 
mass. This is a lower limit since there may be a significant mass 
of H2 too cold to be seen in emission. The H2 kinetic energy is 
more than 3 times the amount of thermal energy of the postshock 
hot plasma. This implies that most of the collision energy is not 
dissipated in the hot plasma, but a substantial amount is carried 
by the warm molecular gas. 

Radio synchrotron observations indicate that the postshock 
medium is magnetized. Assuming the equipartition of energy 
between cosmic rays and magnetic energy, the mean magnetic 
field strength in the SQ shocked region is =s 10yuG (Xuet al., 
2003). The magnetic energy and thereby the cosmic ray energy 
contained in the volume is §- V${ = 2.7 x 10 54 erg. This 
energy is much smaller than the bulk kinetic energy of the H2 
gas. 

To conclude, the dominant postshock energy reservoir is the 
bulk kinetic energy of the molecular gas. We thus propose that 
the observed H2 emission is powered by the dissipation of this 
kinetic energy. The ratio between the kinetic energy of the warm 
H2 gas in the center of mass frame and the H2 luminosity is: 

o/1.5x 10 7 L \ 
fdiss-10 8 — -j yr. (6) 

The H2 luminosity of 1 .5 x 10 7 L is summed over the S(0)-S(5) 
lines. This dissipation timescale is of the order of the cooling 
time of the hot X-ray emitting plasma . Therefore, the dissi- 
pation of the kinetic energy of the SQ collision can power the 
H2 emission on a timescale more than one order of magnitude 
longer than the collision age (f co u ~ 5 x 10 6 yr). The timescale of 
energy dissipation will be larger if the total amount of molecular 
gas is larger than the warm H2 mass derived from Spitzer obser- 
vations. On the other hand, it will be smaller by a factor 2 if we 
consider cooling through the emission of the Hn gas. 



The line emission of the ionized gas could also contribute 
to the dissipation of energy. We estimate the total luminosity of 
the postshock Hn gas from the Ha luminosity (Heckman et al., 
1989): 

13.6 eV 

LHII -0.45x£ H / Hff ' (?) 

where En a is the energy of an Ha photon. Therefore, the lumi- 
nosity of the postshock Hn gas is Lhu = 2.1 x 10 7 L within the 
aperture 3\. This represents 16 x L Ha and is comparable to the 
H2 luminosity. 

3. Why is H 2 present in the postshock gas? 

Observations show that the pre-collision medium is inhomoge- 
neous. The gas in the SQ group halo is observed to be structured 
with Hi clouds and possibly H2 clouds embedded in a teneous, 
hot, intercloud plasma (section 2.1). On the intruder side, the 
galaxy interstellar medium is also expected to be multiphase. 
We consider a large-scale collision between two flows of multi- 
phase dusty gas (Fig. 1). In this section we quantify H2 formation 
within this context. The details of the microphysics calculations 
of the dust evolution, gas cooling and H2 formation are given in 
Appendices A, B and C. 

3. 1 . Collision between two dusty multiphase gas flows 

This section introduces a schematic description of the collision 
that allows us to quantify the H2 formation timescale. The col- 
lision drives two high-velocity large-scale shocks that propa- 
gate through the teneous volume filling plasma, both into the 
tidal arm and into the intruder. The rise in the intercloud pres- 
sure drives slower shocks into the clouds (e.g. McKee & Cowie, 
1975). The velocity of the transmitted shock into the clouds, V c , 

is of order V c - yjrthjn c Vh , where Vh is the shock velocity in 
the hot intercloud medium, and n c and «h are the densities of the 
cloud and intercloud medium 3 , respectively. Each cloud density 
corresponds to a shock velocity. 

The galaxy collision generates a range of shock velocities 
and postshock gas temperatures. The state of the postshock gas 
is related to the preshock gas density. Schematically, low density 
(n H < 10- 2 cnr 3 ) gas is shocked at high speed (> 600 km s 1 ) 
and accounts for the X-ray emission. The postshock plasma did 
not have time to cool down significantly and form molecular gas 
since the collision was initiated. Hi gas («h S 3 x 10~ 2 ) is heated 
to lower (< 10 6 K) temperatures. This gas has time to cool. Since 
H2 forms on dust grains, dust survival is an essential element 
of our scenario (see Appendix A). The gas pressure is too high 
for warm neutral Hi to be in thermal equilibrium. Therefore, the 
clouds cool to the temperatures of the cold neutral medium and 
become molecular. 

The time dependence of both the gas temperature and the 
dust-to-gas mass ratio, starting from gas at an initial postshock 
temperature T v% , has been calculated. We assume a galactic dust- 
to-gas ratio, a solar metallicity and equilibrium ionization. In 
the calculations, the thermal gas pressure is constant. As the 
gas cools, it condenses. Below ~ 10 4 K, the Hn gas recombines, 
cools through optical line emission, and becomes molecular. The 
gas and chemistry network code (Appendix C) follows the time 
evolution of the gas that recombines. The effect of the magnetic 
field on the compression of the gas has been neglected in the 
calculations and we assume a constant thermal gas pressure. 

3 More detailed expressions are given in Klein et al. (1994) 
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Fig. 1. A schematic picture of a galactic wide shock in an inhomogeneous medium. Left: A three data sets image of the Stephan's 
Quintet system (visible red light (blue), Ha (green) from the Calar Alto Observatory in Spain, and 8-micron infrared light (red) from 
Spitzer's InfraRed Array Camera (IRAQ. The white box represents the aperture J{ where Spitzer spectroscopic observations were 
performed. Right: Schematic view of the proposed scenario for the H2 formation in the SQ galaxy-wide collision. In the hot/tenuous 
medium, the shock is propagating fast (V/, ~ 600 km s~'), whereas in the clouds a lower velocity shock (V c < Vh) is driven. The 
cloud compression triggers gas cooling and H2 formation in the postshock clouds. Thermal instability and dynamical interaction 
with the hot plasma lead to molecular gas fragmentation (see § 3.3). 



Under the usual flux-freezing ideas, the magnetic field 
strength should increase with volume density and limit the 
compression factor. However, observations show that the mag- 
netic field in the cold neutral medium (CNM) is compara- 
ble to that measured in lower density, warm ISM components 
(Heiles & Troland, 2005). 3D MHD numerical simulations of 
the supernovae-driven turbulence in the ISM also show that gas 
compression can occur over many orders of magnitude, with- 
out increasing the magnetic pressure (e.g. Mac Low et al., 2005; 
de Avillez & Breitschwerdt, 2005b). The reason is that the lower 
density («h 55 10 3 cirT 3 ) gas is magnetically sub-critical and 
mainly flows along the magnetic field lines without compress- 
ing the magnetic flux. The magnetic field intensity increases 
within gravitationally bound structures at higher gas densities 
(«h <; 10 3 cm 3 , (e.g. Hennebelle et al., 2008). This result is in 
agreement with observations of the magnetic field intensity in 
molecular clouds by Crutcher (1999) and Crutcher et al. (2003). 

3.2. Gas cooling, dust destruction and Hi formation 
timescales 

The physical state of the post-shock gas depends on the shock 
age, the gas cooling, the dust destruction, and the H2 formation 
timescales. These relevant timescales are plotted as a function of 
the postshock temperature in Fig. 2. Given a shock velocity Vh in 
the low-density plasma, each value of the postshock temperature 
T ps corresponds to a cloud preshock density « c p : 

r ps = 2.4xl0 5 f * ,)f°- 1Cm ' 3 )f Vh .fK. (8) 
ps \0.02cm- 3 /\ « c>p /UoOkms-'/ V ; 



Dust destruction is included in the calculation of the H2 for- 
mation timescale. The H2 formation timescale is defined as the 
time when the H2 fractional abundance reaches 90% of its final 
(end of cooling) value, including the gas cooling time from the 
postshock temperature. The dust destruction timescale is defined 
as the time when 90% of the dust mass is destroyed and returned 
to the gas. At T < 10 6 K, this timescale rises steeply because the 
dust sputtering rate drops (see Appendix A for details). The dust 
destruction timescale increases towards higher temperatures be- 
cause, for a fixed pressure, the density decreases. 

The present state of the multiphase postshock gas is set 
by the relative values of the three timescales at the SQ colli- 
sion age, indicated with the horizontal dashed line and marked 
with an arrow. The ordering of the timescales depends on the 
preshock density and postshock gas temperature. By compar- 
ing these timescales one can define three gas phases, marked by 
blue, grey and red thick lines in Fig. 2. 



1. The molecular phase : For r ps < 8 x 10 5 K, the dust de- 
struction timescale is larger than the collision age and the 
H2 formation timescale is lower. This gas keeps a significant 
fraction of its initial dust content and becomes H2 (blue thick 
line). 

2. Atomic Hi and ionized Hn gas phase : For intermediate tem- 
peratures (8 x 10 5 < T ps < 3 x 10 6 K), the gas has time to cool 
down but looses its dust content. The H2 formation timescale 
is greater than the collision age. This phase is indicated with 
the grey thick line. 
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Fig. 2. The multiphasic postshock medium. The dust lifetime, 
H2 formation and gas cooling (solid line) time scales in the SQ 
postshock gas are plotted as a function of the postshock tempera- 
ture and preshock density (top axis). The state of the multiphase 
postshock gas is set by the relative values of these timescales at 
the SQ collision age (~ 5 x 10 6 yr, indicated with the horizon- 
tal dashed line and marked with an arrow). By comparing these 
timescales one can define three gas phases, marked by blue, grey 
and red thick lines. In particular, the comparison of the IT for- 
mation time scale with the SQ collision age show where H2 
molecules can form (blue thick line). 

3. X-ray emitting plasma : At high temperatures (T ps > 3 x 
10 6 K), the postshock gas is dust-free and hot. This is the 
X-ray emitting plasma indicated by the red thick line. 

In Fig. 2 it is assumed that, during the gas cooling, 
the molecular gas is in pressure equilibrium with the hot, 
volume filling gas. The gas pressure is set to the mea- 
sured average thermal pressure of the hot plasma (HIM). 
Theoretical studies and numerical simulations show that the 
probability distribution function of the density in a supersoni- 
cally turbulent isothermal gas is lognormal, with a dispersion 
which increases with the Mach number (Vazquez-Semadeni, 
1994; Padoanetal., 1997; Passot & Vazquez-Semadeni, 1998; 
Mac Low et al., 2005). Based on these studies, the SQ post- 
shock pressure in the turbulent cloud phase is likely to be 
widely distributed around the postshock pressure in the HIM 
(2 x 10 5 K cm -3 ). The H2 formation timescale would be larger 
in the low pressure regions and shorter in high pressure regions. 

To quantify this statement, we have run a grid of isobaric 
cooling models for different gas thermal pressures. Fig. 3 shows 
the gas cooling time and the H2 formation timescale as a func- 
tion of the gas thermal pressure. A Solar Neighborhood value 
of the dust-to-gas ratio is assumed. The H2 formation time scale 
roughly scales as the inverse of the gas density at 10 4 K and can 
be fitted by the following power-law: 



/2x 10 s [KcirT 3 p aS 
tn 2 [yr] - 7xl0 5 f dust 1 



(9) 



where P± is the thermal gas pressure and /dust the dust mass 
fraction remaining in the gas (In Fig. 3, /dust = 1). 

Fig. 3 shows that for a wide range of gas thermal pressures 
(> 3 x 10 4 K cm" 3 , i.e. P th (CNM) > ^P th (HIM)), H 2 can form 
within less than the collision age (5 x 10 6 yrs). Numerical sim- 
ulations of supernovae-driven turbulence in the ISM show that 
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Fig. 3. Cooling and H2 formation timescales for ionized gas at 
10 4 K as a function of the gas density and pressure (top axis). 
The grey horizontal bar indicates the SQ collision age, and the 
blue vertical bar shows the gas thermal pressure at which the H2 
formation timescale equals the collision age. The vertical green 
bar indicates the thermal pressure of the hot plasma. A power- 
law fit to the H2 formation timescale is indicated (Eq. 9). The H2 
formation timescale is shorter than the collision age for pressures 
larger than 1 /10 of the hot plasma pressure. 



this condition applies to a dominant fraction of the gas mass 
(see Fig. 10 of Mac Low et al., 2005). Numerical studies by 
Glover & Mac Low (2007) show that the pressure variations in- 
duced by the turbulence globally reduces the time required to 
form large quantities of H2 within the Galactic CNM. 

3.3. Cloud crushing and fragmentation timescales 

In § 3.2, the microphysics of the gas has been quantified, ig- 
noring the fragmentation of clouds by thermal, Rayleigh-Taylor 
and Kelvin-Helmholtz instabilities. In this section, we introduce 
the dynamical timescales associated with the fragmentation of 
clouds and compare them with the H2 formation timescale. 

3.3.1. Fragmentation 

When a cloud is overtaken by a strong shock wave, a flow of 
background gas establishes around it, and the cloud interacts 
with the low-density gas which is moving relative to it (see 
Fig. 1). The dynamical evolution of such a cloud has been de- 
scribed in many papers (e.g. Klein et al., 1994; Mac Low et al., 
1994). Numerical simulations show that this interaction trig- 
gers clouds fragmentation. The shock compresses the cloud on 
a crushing time ? crush = R c /V c , where V c is the shock velocity 
in the cloud and R c its characteristic preshock size. The shocked 
cloud is subject to both Rayleigh-Taylor and Kelvin-Helmholtz 
instabilities. 

Gas cooling introduces an additional fragmentation mecha- 
nism because the cooling gas is thermally unstable for temper- 
atures smaller than ~ 10 6 K. The thermal instability generates 
complex inhomogeneous structures, dense regions that cool and 
low density voids (Sutherland et al., 2003; Audit & Hennebelle, 
2005). This fragmentation occurs on the gas cooling timescale. 
In most numerical simulations, the thermal instability is ignored 
because simplifying assumptions are made on the equation of 
state of the gas (adiabatic, isothermal or polytropic evolution). 
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Fig. 4. The gas cooling and H2 formation timescales are com- 
pared to the cloud crushing time and SQ collision age as a func- 
tion of the cloud parameters (size and preshock density). The 
solid black line indicates where the gas cooling time f coo i equals 
the crushing time f clus h. The black dashed line show where the 
cloud crushing timescale equals 5 x 10 6 yr (i.e. the SQ collision 
age). The red dashed line defines the location where the H2 for- 
mation timescale equals the crushing timescale. The red area 
indicates where the H2 formation timescale is greater than the 
collision age, whereas the blue hatched zone show the parame- 
ter space where H2 can form. The yellow zone is the location of 
gravitationally bound, unstable clouds that have a free-fall time 
shorter than the collision age. 



In Fig. 4, the Ht formation and gas cooling timescales are 
compared to the crushing time and collision age, as a function 
of the initial density and size of the cloud. The pressure of the 
clouds is 2 x 10 5 K cm 3 . The red zone indicates where the H2 
formation timescale is larger than the SQ collision age. In this 
domain, clouds do not have time to form H2. The blue hatched 
zone defines the range of sizes and densities of the clouds which 
become molecular within 5 x 10 6 yr. In this zone, the crushing 
time is larger than the cooling time. Gas fragments by thermal 
instability as the shock moves into the cloud. H2 forms before 
the shock has crossed the cloud and therefore before Rayleigh- 
Taylor and Kelvin-Helmholtz instabilities start to develop. The 
black dashed line is where the cloud crushing timescale equals 
the SQ collision age (5 x 10 6 yr). To the right of this line, the 
clouds have not been completely overtaken and compressed by 
the transmitted shock. 

Preshock Giant Molecular Clouds (hereafter GMCs) of size 
> 30 pc and mean density > 3 x 10 2 cm 3 , like those observed in 
spiral galaxies, fall within the yellow zone in Fig. 4. The bottom 
boundary corresponds to the density at which the free-fall time 
equals the age of the SQ collision. For densities higher than this 
limit, the free-fall time is shorter than the SQ collision age. The 
left boundary of this area is the maximum cloud size set by the 
Bonnor-Ebert criterion for a pressure of2x 10 5 K cm . Clouds 
with larger sizes are gravitationally unstable. Such clouds are 
expected to form stars. 

A natural consequence of pre-existing GMCs would be that 
the shock would trigger their rapid collapse and star formation 
should rapidly follow (e.g. Jog & Solomon, 1992). The absence 
(or weakness) of these tracers in the region of the main shock 
(Xu et al., 2003) suggests that fragments of preshock GMCs dis- 
rupted by star formation are not a major source of the observed 



postshock H2 gas. However, it is possible that some of the gas 
in the intruder did contain pre-existing molecular material. The 
star formation regions seen to the North and South of the ridge 
are consistent with this hypothesis. 

3.3.2. Survival of the fragments 

The breakup of clouds can lead to mixing of cloud gas with 
the hot background plasma (Nakamura et al., 2006) and thereby 
raises the question whether the cloud fragments survive long 
enough for H2 to form. Fragments evaporate into the back- 
ground gas when their size is such that the heat conduction 
rate overtakes the cooling rate. The survival of the fragments 
thus depends on their ability to cool (Mellema et al., 2002; 
Fragile et al., 2004). Gas condensation and H2 formation are 
both expected to have a stabilizing effect on the cloud fragments, 
because they increase the gas cooling power (Le Bourlot et al., 
1999). 

So far, no simulations include the full range of densities and 
spatial scales and the detailed micro-physics which are involved 
in the cloud fragmentation and eventual mixing with the back- 
ground gas. Slavin (2006) reports an analytical estimate of the 
evaporation time scale of tiny H 1 clouds embedded in warm and 
hot gas of several million years. This timescale is long enough 
for H2 to form within the gas fragments lifetime. 



4. H 2 excitation 

This section addresses the question of H2 excitation. The models 
(§ 4.1) and the results (§ 4.2) are presented. 

4.1. Modeling the dissipation of turbulent kinetic energy by 
MHD shocks 

Given that the dominant postshock energy reservoir is the bulk 
kinetic energy of the H 2 gas (§ 2.4), the H2 emission is assumed 
to be powered by the dissipation of turbulent kinetic energy into 
molecular gas. Radio synchrotron observations indicate that the 
postshock medium is magnetized (§ 2.4). In this context, the dis- 
sipation of mechanical energy within molecular gas is modeled 
by non-dissociative MHD shocks. In the absence of magnetic 
field, shocks are rapidly dissociative and much of the cooling 
of the shocked gas occurs through atomic, rather than H2, line 
emission. We quantify the required gas densities and shock ve- 
locities to account for the H2 excitation diagram. 

The H2 emission induced by low velocity shocks is cal- 
culated using an updated version of the shock model of 
Flower & Pineau des Forets (2003). In these models, we assume 
a standard value for the cosmic ray ionization rate of ( = 
5 x 10~ 17 s _1 . In its initial state, the gas is molecular and cold 
(~ 10 K), with molecular abundances resulting from the out- 
put of our model for the isobaric cooling (§ C.2 and Fig. C.2). 
A grid of shock models has been computed for shock veloci- 
ties from 3 to 40 km s~' with steps of 1 km s _1 , two values of 
the initial H2 ortho to para ratio (3 and 10~ 2 ), and 3 different 
preshock densities («h = 10 2 , 10 3 , and 10 4 cm 4 ). We adopt 
the scaling of the initial magnetic strength with the gas den- 
sity, B(/jG) - b V«h [cm^ 3 ] (Crutcher, 1999; Hennebelle et al., 
2008). In our models, b = 1. 
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Fig. 5. H2 excitation diagram fit by shock models. The data (Appleton et al., 2006) are the black crosses with their l-cr error bars. 
The blue dashed line show the modeled excitation associated with the H2 formation during the isobaric cooling of the postshock 
gas. The pressure is 2 x 10 5 K cm 4 . The model column densities correspond to a mass flow of recombining gas of 100 M Q yr . 
The green and orange dashed lines show the contribution of 5 and 20 km s shocks to the best fit model (red line), multiplied 
by the weights determined when fitting the line fluxes, respectively. Shocks are driven into molecular gas of preshock density 
«h = 10 4 cm" 3 . The red line shows the linear combination of the two magnetic shocks (with velocities of 5 and 20 km s _1 ). 



4.2. Results: contribution of MHD shocks to the H 2 excitation 

Results are presented in terms of an H2 excitation diagram, 
where the logarithm of the column densities N v j[cm~ 2 ], divided 
by the statistical weight g Y j, are plotted versus the energy of the 
corresponding level. The excitation diagrams are calculated for 
a gas temperature of 50 K, since the molecular gas colder than 
50 K does not contribute to the H2 emission. Fig. 5 shows the 
H2 excitation diagram for both observations and model. The SQ 
data from Appleton et al. (2006) is indicated with lcr error bars. 
The H2 line fluxes have been converted to mean column densities 
within our area M. 

The observed rotational diagram is curved, suggesting that 
the H2-emitting gas is a mixture of components at different tem- 
peratures. The data cannot be fitted with one single shock. At 
least two different shock velocities are needed. A least-squares 
fit of the observed H2 line fluxes to a linear combination of 2 
MHD shocks at different velocities V s i and V S 2 has been per- 
formed. For a preshock density of 10 4 cm -3 , the best fit is ob- 
tained for a combination of 5 and 20 km s _1 shocks with an ini- 
tial value of the ortho to para ratio of 3 (Fig. 5 and Table 2). Note 
that no satisfactory fit could be obtained with the low value of the 



initial H2 ortho-to-para ratio. The contributions of the 5 (green 
dashed line) and 20 km s _1 (orange dashed line) shocks to the to- 
tal column densities are also plotted. Qualitatively, the 5 km s _1 
component dominates the contribution to the column densities 
of the upper levels of S(0) and S(l) lines. The 20 km s com- 
ponent has a major contribution for the upper levels of S(3) to 
S(5). For S(2), the 2 contributions are comparable. The red line 
shows the weighted sum of the two shock components (best fit). 
The warm H2 pressure p V 2 in the two shocks is 4.5 x 10 7 K cirT 3 
(respectively 5.6x1 8 Kern -3 ) for the 5 km s _1 (resp. 20 km s _1 ) 
shock. 

Our grid of models shows that this solution is not unique. 
If one decreases the density to 10 3 cirT 3 , we find that one can 
fit the observations with a combination of MHD shocks (at 
9 and 35 km s _1 ). If one decreases the density to 10 2 cm -3 , 
the rotational H2 excitation cannot be fitted satisfactorily. At 
such low densities MHD shocks fail to reproduce the S(3) and 
S(5) lines because the critical densities for rotational H2 exci- 
tation increases steeply with the J rotational quantum number 
(Le Bourlot et al., 1999). Both warm and dense (i.e. high pres- 
sure) gas is needed to account for emission from the higher J 
levels. 
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The data fit gives an estimate of the shock velocities required 
to account for the H2 excitation. The velocities are remarkably 
small with respect to the SQ collision velocity. The fact that en- 
ergy dissipation occurs over this low velocities range is an es- 
sential key to account for the importance of fT cooling. High 
velocity shocks would dissociate H2 molecules. 

5. The cycling of gas across ISM phases 

An evolutionary picture of the inhomogeneous postshock gas 
emerges from our interpretation of the data. The powerful H2 
emission is part of a broader astrophysical context, including 
optical and X-ray emission which can be understood in terms of 
mass exchange between gas phases. A schematic cartoon of the 
evolutionary picture of the postshock gas is presented in Fig. 6. 
It illustrates the cycle of mass between the ISM phases that we 
detail here. This view of the SQ postshock gas introduces a phys- 
ical framework that may apply to IT luminous galaxies in gen- 
eral. 

5.1. Mass flows 

The H2 gas luminosity is proportional to the mass of gas that 
cools per unit time, the so-called mass flow. We compute the 
mass flow associated with one given shock dividing the post- 
shock column density down to a gas temperature of 50 K by the 
cooling time needed to reach this temperature. Gas cooler than 
50 K does not contribute to the H2 emission. 

The two-components shock model has been used to deduce 
the mass flows needed to fit the data. The cooling timescale 
(down to 50 K) is set by the lowest shock velocity component. 
For a shock velocity of 5 km s , the cooling time is ~ 2 x 10 4 yr. 
This gives a total mass flow required to fit the IT emission of 
M3 = 5.7 x 1O 3 M yr -1 . For comparison, for the 20 km s _1 
shock, the cooling time down to 50 K is ~ 4 x 10 3 yr and the 
associated mass flow is M = 270 M Q yr -1 . 

These mass flows are compared to the mass flow of recom- 
bining gas that can be estimated from the H a luminosity Lu a ■ 

L Ha = 0.45 x E Ho x <t> Rec , (10) 

where En a is the energy of an H a photon and <t>R ec is the number 
of hydrogen atoms that recombine per second (Heckman et al., 
1989). For our reference area M, we find a mass flow of recom- 
bining gas Mu„ = '«h 'I'Rec — 100 M Q yr -1 . If the ionizing rate 
of warm gas is smaller than <J>Rec, most of the recombining gas 
cools. In this case <I>R ec provides the estimate of the mass flow 
needed to compute the cooling gas H2 luminosity. Or cc is an up- 
per limit if we suppose that some of the H ff emission arises from 
collisional excitation of neutral hydrogen or if some of the warm 
gas is re-ionized before it has time to cool. 

The H2 emission integrated over the isobaric cooling se- 
quence is calculated. The blue dashed curve of Fig. 5 corre- 
sponds to the model results for a mass flow of M = 100 M yr -1 
and a pressure Psq = 2 x 10 5 K cm 4 . The model results do 
not depend much on the pressure. The model reproduces well 
the shape of the observed excitation diagram, but the emission 
is a factor ~ 25 below the observations. The model also fails 
to reproduce the [Oi] 6 300 A line emission by a factor of 7. 
Therefore, the cooling of the warm Hi gas produced by H11 re- 
combination cannot reproduce the observed H2 emission nor the 
[Oi] 6 300 A line emission. A larger amount of energy needs 
to be dissipated within the molecular gas to account for the H2 
emission. This is why the fit total mass flow associated with the 



Mass and energy cycle 




Fig. 6. Schematic view at the evolutionary cycle of the gas pro- 
posed in our interpretation of optical and H2 observations of 
Stephan's Quintet. Arrows represent the mass flows between the 
Hn, warm Hi, warm and cold H2 gas components. They are num- 
bered for clarity (see text). The dynamical interaction between 
gas phases drives the cycle. The values of the mass flows and 
associated timescales are derived from the Ha, [Oi 6 300 A], 
H2 observed luminosities and model calculations (cf. Table 2). 
Heating of the cold H2 gas (red dashed arrows) is necessary to 
account for the increasing mass flow from the ionized gas to cold 
H2 phases. 



2-component MHD shocks is much larger than the mass flow of 
recombining gas. 

5.2. The mass cycle 

This section describes the mass cycle of Fig. 6. Black and red ar- 
rows represent the mass flows between the Hn, warm Hi, warm 
and cold H2 gas components of the postshock gas. The large 
red arrow to the left symbolizes the relative motion between 
the warm and cold gas and the surrounding plasma. Each of the 
black arrows is labeled with its main associated process: gas re- 
combination and ionization (double arrow number 1), H2 forma- 
tion (2) and H2 cooling (3). The values of the mass flows and 
the associated timescales are derived from observations and our 
model calculations (§5.1 and Table 2). 

A continuous cycle through gas components is excluded by 
the increasing mass flow needed to account for the Ha, Oi, and 
H2 luminosities. Heating of the cold H2 gas towards warmer gas 
states (red arrows) must occur. The postshock molecular cloud 
fragments are likely to experience a distribution of shock veloc- 
ities. Arrow number 4 represents the low velocity MHD shock 
excitation of H2 gas described in § 4. 1 . A cyclic mass exchange 
between the cold molecular gas and the warm gas phases is also 
supported by the high value of the H2 ortho to para ratio (§ 4.1). 
More energetic shocks may dissociate the molecular gas (arrows 
number 5). They are necessary to account for the low O 1 lu- 
minosity. Even more energetic shocks can ionize the molecu- 
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Table 2. Overview of the emission of the SQ postshock gas, from observations and model predictions. All line fluxes are scaled to 
the aperture J{ — 1 1 .3 x 4.7 arcsec 2 . For isobaric cooling calculations, the pressure is set to the SQ postshock gas value. For MHD 
shock models, both contributions of the 2 shock velocities to the emission are indicated, as well as the pressure for the warm H2 
phase in the shock. These models are the same as those used in Fig. 5 



" MHD Shock velocity. 

Cooling time computed down to 50 K. 
c Sum of the H 2 S(0) to S(5) rotational lines (from Appleton et al. (2006)). 
d From optical observations by Xu et al. (2003). 

e Mass flow derived from H a observations, ignoring ionization. It represents the isobaric cooling of the recombining Hn gas. 
' Cooling mass flow derived from Oi observations and model calculations at the ambient SQ pressure. This mass flow does not depend much on 
the pressure since the critical density of the Oi line is high (2 x 10 6 cm -3 ). 



lar gas (arrow number 6). This would bring cold H2 directly 
into the H 11 component. Such shocks have been proposed as 
the ionization source of the clouds in the Magellanic Stream by 
Bland-Hawthorn et al. (2007) based on hydrodynamical simula- 
tions of their interaction with the Galactic halo plasma. 

Turbulent mixing of hot, warm and cold gas are an alterna- 
tive mass input of H n gas. Turbulent mixing layers have been de- 
scribed by Begelman & Fabian (1990) and Slavin et al. (1993). 
They result from the shredding of the cloud gas into fragments 
that become too small to survive evaporation due to heat con- 
duction. It is represented by the thin blue and green arrows to 
the left of our cartoon. Turbulent mixing produces intermediate 
temperature gas that is thermally unstable. This gas cools back 
to produce H 11 gas that enters the cycle (thin green arrow). It 
is relevant for our scenario to note that cold gas in mixing lay- 
ers preserves its dust content. It is only heated to a few 10 s K, 
well below temperatures for which thermal sputtering becomes 
effective. Further, metals from the dust-free hot plasma brought 
to cool are expected to accrete on dust when gas cools and con- 
denses. Far-UV Ovi line observations are needed to quantify the 
mass flow rate through this path of the cycle and thereby the rel- 
ative importance of photoionization (arrow 1), ionizing shocks 
(arrow 6), and turbulent mixing with hot gas in sustaining the 
observed rate (^ 100 M yr -1 ) of recombining H 11 gas. 

Within this dynamical picture of the postshock gas, cold 
molecular gas is not necessarily a mass sink. Molecular gas 
colder than 50 K does not contribute to the H2 emission. The 
mass accumulated in this cold H2 gas depends on the ratio be- 
tween the timescale of mechanical energy dissipation (on which 
gas is excited by shocks through one of the red arrows 4, 5 and 
6) and the cooling timescale along the black arrow 3. 

6. Why H 2 is a dominant coolant of the postshock 
gas? 

This section sketches the energetics of the molecular gas in rela- 
tion to its interaction with the background plasma. An efficient 
transfer of the bulk kinetic energy of the gas to turbulent motions 
of much lower velocities within molecular clouds is required to 
make H2 a dominant coolant of the postshock gas. We consider 



that the bulk kinetic energy of the clouds drives the gas cycling 
through the warm/cold H n, H 1 and H2 states of multiphase frag- 
mentary clouds and feeds the cloud turbulence. This hypothesis 
is supported by the fact that the mechanical energy is the main 
energy reservoir (§ 2.4). The energetics of the gas are discussed 
within this framework. 

6. 1. Cloud dynamics and energy transfer 

In the collision, the clouds do not share the same dynamics as 
the lower density intercloud gas. The clouds keep their preshock 
momentum and move at high velocity with respect to the back- 
ground plasma. The kinetic energy recorded in the 900 km s _1 
H2 line width is dominated by the bulk motion of the H2 clouds 
in the plasma rest frame. The turbulent velocities required within 
clouds to excite Fb have a much smaller amplitude (relative frag- 
ment velocities up to 20 kms 1 ). The cloud turbulence needs to 
be continuously sustained to balance energy dissipation. An ef- 
ficient transfer of the bulk kinetic energy to low velocity turbu- 
lence of the molecular gas is required to make H2 a dominant 
coolant of the postshock gas. Dissipation within either the hot 
plasma or warm atomic gas, leading to radiation at X-ray or opti- 
cal wavelengths that is not efficiently absorbed by the molecular 
gas, has to be minor. 

We are far from being able to describe how the energy trans- 
fer occurs. We just make qualitative statements that would need 
to be quantified with numerical simulations in future studies. 
Cloud motions are subject to drag from the background flow 
(Murray & Lin, 2004). This drag corresponds to an exchange 
of momentum and energy with the background plasma which 
drives turbulence into the clouds. 

The fragmented, multiphase structure of the clouds is likely 
to be a main key of the energy transfer. The dynamical interac- 
tion with the background plasma flow generates relative veloci- 
ties between cloud fragments because their acceleration depends 
on their mass and density. The relative motions between cloud 
fragments can lead to collisions, which results in the dissipation 
of kinetic energy (§ 4.1). This dissipation will occur preferen- 
tially in the molecular gas because it is the coldest component 
with the lowest sound speed. The magnetic field is likely to be 
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important because it weaves the cloud H2 fragments to the lower 
density cloud gas which is more easily entrained by the back- 
ground plasma flow. The mass cycle between cloud gas states 
contributes to the energy transfer in two ways. (1) Gas cooling 
transfers the turbulent energy of the warm H 11 and H 1 to the H2 
gas. (2) The thermal instability which occurs when the gas cools 
both from 10 6 K and from 10 4 K transfers a significant fraction 
of the gas thermal energy to turbulent kinetic energy. Numerical 
simulations (Sutherland et al., 2003; Audit & Hennebelle, 2005) 
show that the thermal instability inhibits energy-loss through ra- 
diation from the warmer gas and feeds turbulence in the cooler 
gas. The turbulent velocities generated by the thermal instabil- 
ity are found to be commensurate with the sound speed in the 
warmer gas and are thus supersonic for the colder gas. 

6.2. Gas energetics 

The cloud turbulent energy and the dissipated power per H2 
molecule may be written as 

3 9 

E T = -m Hl o-^ T (11) 

P T = ?L = 0.24(-^) 2 (-^-\ 1 l3>, (12) 
t T UOkms" 1 / \10 5 yr/ M 

where toh, is the mass of the H2 molecule and crj the cloud 
turbulent velocity dispersion. The second equation introduces an 
effective dissipation timescale, tj which is often written as the 
ratio between the cloud size, L, and <tt- The reference value of 
<Tt used in the second equation is commensurate with the shock 
velocities inferred from the fitting of the H2 excitation diagram 
in §4.1. 

If dissipation of turbulent energy powers the H2 emission, 
the energy dissipation rate must be at least equal to the power 
radiated by H2. Using the H2 emission and the warm H2 mass 
in Table 1, we estimate the luminosity to mass ratio of the warm 
H2 gas in SQ to be 0.5 L /M o . The dissipated power Pj may 
be somewhat larger if it powers the dissociative and ionizing 
shocks represented by the arrows 5 and 6 in Fig. 6. This cal- 
culation implies a dissipation timescale of ~ 5 x 10 4 yr which 
is not much higher than the cooling time of the lower excitation 
warm H2 gas which accounts for the S(0) and S(l) line emission 
(§ 5.1). Cold H2 gas would thus be heated by turbulent dissi- 
pation on timescales comparable to the gas cooling time. If this 
is right, the cold molecular gas mass should not be much larger 
than the warm H2 mass inferred from the Spitzer observations. 
The turbulence dissipation timescale is more than three orders of 
magnitude smaller than the ratio between the H2 gas bulk kinetic 
energy and H2 luminosity (section 2.1). Thereby, the dissipation 
of the collision energy must involve a large (> 1 000) number 
of cycles where H2 gas fragments are accelerated and dissipate 
their kinetic energy in shocks. 

7. Future tests of the model and open questions 

In this section we provide a few observational predictions of the 
model related to questions that remain open. For each point, we 
discuss which observations may test these predictions. 

7. 1 . Dust Emission 

Dust is a key element in our scenario. It must be present in the 
postshock clouds for H2 to form. Xu et al. (1999, 2003) present 



results from infrared imaging observations carried out with the 
Infrared Space Observatory (ISO). Far-infrared (FIR) emission 
was detected in the shock region and Xu et al. (2003) proposed 
that it arises from large dust grains in the hot X-ray emitting 
plasma. Clearly, the Spitzer detection of molecular gas opens a 
new perpective on the origin of this emission. In our model we 
expect that some of the thermal dust emission comes from the 
molecular gas. Outside local UV enhancement associated with 
star formation seen in GALEX observations (Xu et al., 2005), 
dust emission should appear spatially correlated with the H2 
emission. Given that the value of the radiation field in the shock 
is low (G ~ 1 - 2 Go § 2.3), we expect rather faint dust emission 
in the SQ ridge. This accounts for the weakness of the PAH fea- 
tures in the Spitzer spectra. Deep Spitzer imaging observations 
have been undertaken and will allow us to test this prediction. 
Moreover, observations with the Herschel Space Telescope will 
have the required angular resolution and sensitivity to detect the 
long wavelength emission of the dust in the SQ shock. 



7.2. Mass Budget 

The Spitzer H2 observations only trace the warm molecular gas 
(T > 150 K). This warm gas could be physically connected to 
a more massive reservoir of cold H2 gas. This is a major open 
question about the physical state of the molecular gas, as well 
as about the mass and energy budget. In a forthcoming paper, 
we will report deep IRAM CO observations that set a low up- 
per limit on the CO(1-0) and (2-1) line emission from the SQ 
shock area. The physical state of the shocked gas is likely to be 
different from that of GMCs in galactic disks which are gravita- 
tionally bound. It is thus not clear that the standard CO emission 
to H2 mass conversion factor applies. Herschel observations of 
the dust emission are needed to trace cold gas independently of 
the CO/H 2 ratio. 



7.3. Shock Energetics 

Spitzer observations of the low-excitation H2 rotational lines 
do not provide a complete budget of the kinetic energy dissi- 
pation. To constrain the present scenario, observations of both 
higher and lower excitation lines are required. Rovibrational H2 
lines would probe dissipation in more energetic shocks than 
those inferred from the H2 pure rotational lines. If the molec- 
ular gas only experiences low-velocity MHD shocks, our model 
predicts a 1-0 S(l) line flux of 3.6 x 10~ 19 W irT 2 within the 
aperture J{ (Table 2). This is likely to be a strict lower limit 
since the dynamical interaction between the gas phases may 
well drive higher shock velocities (§ 3.1). Observations of the 
Oi [63 fim] line would provide constraints on the physical state 
of colder molecular gas. The luminosity ratio of Fb rotational 
lines to Oi depends on the shock velocity and the gas density. 
For the SQ ridge, the observed H2 0-0 S(l) line flux within the 
area J{ is 2 x 10~ 17 W m 2 , and our model predicts a H2 0- 
S(l) to Oi [63 pm] luminosity ratio of 2.5. The Integral Field 
Spectrometer of the PACS instrument onboard the Herschel 
satellite should be able to detect the Oi[63 pm] line in one hour 
of integration. The spectral resolution of these observation will 
be a factor ~ 5 higher than Spitzer. The Oi [63 pm] line pro- 
files will allow to compare the kinematics across the different 
temperatures of the gas. If the emission arises from multiphase 
clouds where warm and cold components are mixed down to 
small scales, their kinematics should be the same. 
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7.4. Star Formation in the postshock gas 

Many galaxy collisions and mergers are observed to trigger 
IR-luminous bursts of star formation. However, the absence or 
weakness of spectroscopic signatures of photoionization (dust or 
ionized gas lines) in the center of the SQ ridge (Xu et al., 2003), 
where bright H2 emission was detected by Spitzer, show that 
star formation does not occur in this region. Within the physi- 
cal framework of dynamical interaction between gas phases de- 
scribed in Sect. 5 and 6, these observations suggest that most 
of the H2 emission arises from molecular fragments that are 
not massive enough or not long-lived enough to collapse. Star 
formation may also be quenched by continuous heating of the 
molecular gas and/or by the magnetic pressure if the dense 
cloud fragments are not magnetically supercritical (§ 3.1). Star- 
forming regions possibly associated with the galaxy collision 
are present away from the center of the shock (Xu et al., 1999; 
Gao & Xu, 2000; Smith & Struck, 2001; Lisenfeld et al., 2002; 
Xu et al., 2003). 

In our model, star formation may arise from the collapse 
of pre-existing GMCs (see § 3.3). An alternative route to star 
formation may be the formation of gravitationally unstable 
postshock fragments. This second possibility has been quanti- 
fied by numerical simulations that have introduced the concept 
of turbulence-regulated star formation (Mac Low & Klessen, 
2004). Krumholz & McKee (2005) present an analytical de- 
scription of the results of simulations, where the star forma- 
tion rate depends on the Mach number of the turbulence within 
clouds and the ratio between their gravitational and turbulent ki- 
netic energy. The star formation rate increases with the Mach 
number, but it is low even for high Mach numbers when the ra- 
tio between binding and turbulent energy is low. In the center 
of the SQ ridge, the star formation may be low because the av- 
erage gas density is low (see Table 1), even if the amount of 
turbulent energy is high. The low gas density and large-scale 
dynamical context (high-velocity dispersion) does not favor the 
formation of gravitationnally bound molecular clouds, such as 
galactic GMCs. 

8. Conclusions 

Spitzer observations led to the surprising detection of H2 lumi- 
nous extragalactic sources whose mid-infrared spectra are dom- 
inated by molecular hydrogen line emission. The absence or 
weakness of dust features or lines of ionized gas suggests the 
presence of large quantities of warm molecular gas with no or 
very little star formation. This work focuses on one of these 
H2 luminous galaxies: the Stephan's Quintet (SQ) galaxy-wide 
shock created by the collision between an intruding galaxy and 
a tidal arm at a relative speed ~ 1000 km s" 1 . SQ observations 
place molecular gas in a new context where one has (1) to ac- 
count for the presence of H2 in a galaxy halo, (2) to characterize 
its physical state and (3) to describe its role as a cooling agent 
of an energetic phase of galaxy interaction. The aim of the pa- 
per was to answer three main questions: (;) Why is there H2 in 
the postshock gas ? (ii) How can we account for the H2 exci- 
tation ? ( Hi) Why is H2 the dominant coolant ? We summarize 
the main results of our work along these three points and present 
perspectives for future studies. 

We have presented and quantified a scenario where H2 forms 
out of shocked gas. H2 formation results from the density struc- 
ture of the preshock gas. The collision velocity, ~ 600 km s" 1 in 
the centre of mass frame of the collision, is the shock velocity 
in the low density («h S5 10~ 2 cm 4 ) volume filling gas. This 



produces a ~ 5 x 10 6 K, dust-free, X-ray emitting plasma. The 
pressure of this background plasma drives slower shocks into 
denser gas within clouds. Gas heated to temperatures smaller 
than ~ 10 6 K keeps its dust and has time to cool within the 
5 x 10 6 yr collision age. Because the postshock pressure is too 
high for the warm neutral medium to be thermally stable, the gas 
cools down to cold neutral medium temperatures, condenses and 
becomes molecular. We show that for a wide range of postshock 
pressures, as well as realistic preshock cloud sizes and densities, 
the H2 formation timescale is shorter than the collision age and 
also than the turbulent mixing timescale of the warm gas with 
the background plasma. 

Observations show that the H2 non-thermal kinetic energy is 
the dominant energy reservoir. It is larger than the thermal en- 
ergy of the plasma. We propose that the H2 emission is powered 
by the dissipation of this kinetic energy. The H2 excitation can be 
reproduced by a combination of low velocity (~ 5 to 20 km s" 1 ) 
MHD shocks within dense («h > 10 3 cm" 3 ) IT gas. In this in- 
terpretation, the pressure of the warm IT gas is much higher than 
the ambient pressure set by the background plasma. Such a pres- 
sure enhancement is required to account for the higher excitation 
H2 S(3) and S(5) rotational lines. For this shock velocity range, 
the warm H2 pressure in the shocks is ~ 0.5 - 5 x 10 8 K cm" 3 , 
2-3 orders of magnitude greater than the ambient pressure 
(2 x 10 5 K cm 4 ). 

The thermal, Rayleigh-Taylor and Kelvin-Helmhotz instabil- 
ities contribute to produce fragmentary postshock clouds where 
H2, Hi and Hn gas are mixed on small scales and embedded 
in the hot plasma. In our scenario, the H2 excitation is related 
to the dynamical interaction between gas phases which drives 
an energy transfer between these phases. An efficient transfer 
of the kinetic energy associated with the bulk displacements of 
the clouds to turbulent motions of much lower velocities within 
molecular gas, is required to make H2 a dominant coolant of the 
postshock gas. 

We present a global view at the postshock gas evolution that 
connects H2, warm Hi and Hn gas emission to mass cycling be- 
tween these gas components. Within this dynamical picture of 
the postshock gas, cold molecular gas is not necessarily a mass 
sink. The fact that there is no star formation in the shocked re- 
gion where H2 is detected suggests that molecular fragments are 
not long-lived enough to collapse and become gravitationally un- 
stable. 

This global description provides a physical framework of 
general relevance for the interpretation of observational signa- 
tures, in particular H2 emission, of mechanical energy dissipa- 
tion in multiphase gas. This study highlights the need to take 
into account H2 formation and the mechanical energy dissipa- 
tion within cold molecular gas in the physics of turbulent mixing 
of gas phases. 

Future observations of the Stephan's Quintet will contribute 
to test the proposed interpretation and better characterize some 
of its aspects. Some of these observations have been briefly dis- 
cussed. Further work is needed to explore the impact of key 
physical parameters (gas pressure, radiation field, mechanical 
energy deposition, thermal to mechanical energy ratio, metal- 
licity) on gas evolution and on energy dissipation for the di- 
versity of relevant galaxy-wide environments and mechanical 
energy sources (gas accretion, galaxy interaction and feedback 
from starburst and active galactic nuclei). The analytical ap- 
proach presented here should be pursued to compute how the 
timescales of the relevant physical processes depend on these 
parameters. The non-linear dynamical interaction between gas 
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phases, in particular the driving of clouds turbulence, can only 
be investigated quantitatively through numerical simulations. 
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Appendix A: Modeling Dust Evolution 

Interstellar dust plays a major role in the evolution of the post- 
shock gas because its survival is required for H2 formation. For 
a constant dust-to-gas mass ratio, the dust dominates the cooling 
efficiency of the gas at high (T > 10 6 K) temperatures (Dwek, 
1987) but for such high temperatures, the timescale for dust de- 
struction is shorter than the gas cooling time (Smith et al., 1996). 
This appendix details how the evolution of the dust-to-gas mass 
ratio in the postshock gas has been calculated. This calculation 
is coupled to that of the gas cooling, detailed in Appendix B. 

The dominant dust destruction process depends on the shock 
velocity. In this study we consider the effect of thermal and iner- 
tial sputtering by grain-ion collisions. This is the dominant mode 
of destruction in a hot (~ 10 6 -10 8 K) plasma resulting from fast 
(100-1 000 km s _1 ) intergalactic shocks (e.g. Draine & Salpeter, 
1979; Dweketal., 1996). 

It is assumed that dust grains have an initial radius equal to 
the effective (mean) radius of a e g = 0.05 //m. This radius has 
been calculated from the MRN dust size distribution (§ B.l). At 
each time step, the fraction of dust remaining in the gas /d us t, is 
computed in parallel to the postshock gas cooling: 



/dust = — I with a 

\fleff; 



a e ff 



I 

Jin 



adt . 



(A.l) 



where a = da/dt. For a shock velocity V s > 300 km s , thermal 
sputtering dominates over inertial sputtering, so we ignore iner- 
tial sputtering. In this case the grain sputtering rate is nearly con- 
stant for T > 3 x 10 6 K and strongly decreases at lower tempera- 
tures (Draine & Salpeter, 1979; Tielens et al., 1994; Dwek et al., 
1996). The rate of decrease in grain size is given by 

da m sn sr~t 

-Z = T L *h. 5>,<r f v> , (A.2) 
dt 2po 1 

where m sp and po are the average mass of the sputtered atoms and 
the density of the grain material, A, the abundance of impacting 
ion i and T, is the sputtering yield 4 of ion i. The quantity (T, v) is 
the average of the impacting ion velocities v times the yield over 
the Maxwellian distribution. The right hand side of this equation 
is computed from the results given in Tielens et al. (1994). 

For V s < 300 km s , we integrate the equation of decelera- 
tion of the dust grains with respect to the gas 

dv g Pna 2 pv 2 g 
dt m 



(A3) 



where v g is the grain velocity with respect to the gas, m and a 
are its mass and radius, p is the gas density, and /3 is the en- 
hancement of the collisional drag in a plasma relative to that of 
a neutral medium. The initial grain velocity is set to 3 /4 of the 
shock velocity (Jones et al., 1994). We take (3 — 1, which max- 
imizes the dust destruction (see Borkowski & Dwek (1995) for 
details). Therefore, our computation of the dust survival is con- 
servative. When inertial sputtering is the dominant destruction 
process, the erosion rate is given by (Tielens et al., 1994) 



A = 2^ v ' BH 2> y '- (A ' 4) 

By integrating these differential equations, the radius of the dust 
grains as a function of their velocity has been deduced (see 
Borkowski & Dwek (1995) for a complete parametrization of 
da/dt in the case of inertial sputtering). 
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Fig. A.l. Dust mass fraction as a function of the initial (post- 
shock) temperature, after isobaric cooling to 10 4 K. The pressure 
is fixed to P ps /Icb = 2.3 x 10 5 [cm~ 3 K]. The dotted line shows 
the result of the model which includes both inertial and thermal 
sputtering. As a comparison, the dashed line only shows the ef- 
fect of thermal sputtering. 

Fig. A.l shows the dust mass fraction remaining in the gas 
after isobaric cooling from the postshock temperature T ps to 
10 4 K as a function of the shock velocity. The dotted line in- 
cludes both thermal and inertial sputtering. For comparison, the 
dashed line only includes dust destruction by thermal sputter- 
ing. Our calculation for thermal sputtering with a single grain 
size (fl e ff = 0.05 pm) is in very good agreement with the study 
by Smith et al. (1996) who used a grain size distribution. For 
T ps > 2 x 10 6 K, i.e. V s > 400 km s" 1 , most of the dust is de- 
stroyed before the gas has cooled down to 10 4 K. The gas which 
is heated to less than ~ 10 6 K keeps a major part of its dust con- 
tent and may therefore form H2. 

Appendix B: Modeling the cooling of hot gas 

This appendix details how we model the cooling of the hot gas. 
The cooling functions used in the calculations are presented in 
§ B.l. Then the results are discussed in § B.2. 

The time dependence of both the gas temperature and the 
dust-to-gas mass ratio has been calculated. We start from gas 
at an initial postshock temperature T ps with a galactic dust-to- 
gas ratio, a solar metallicity and assume equilibrium ionization. 
From a range of postshock temperatures to 10 4 K, the isobaric 
gas cooling is calculated by integrating the energy balance equa- 
tion which gives the rate of decrease of the gas temperature: 



~z &b ~r~ = ~H n e (/d ust A dust + A gas ) , 



(B.l) 



4 The yield is the number of atoms ejected from the grains surface 
per incident particle. 



where p is the mean particle weight (p = 0.6 a.m.u for a fully 
ionized gas), n e is the electronic density, the Boltzman con- 
stant, /j ust the dust-to-gas mass ratio, A dus t and A gas are respec- 
tively the dust and gas cooling efficiencies per unit mass of dust 
and gas, respectively. 

The time-dependent total cooling function is the sum of the 
dust cooling (weighted by the remaining fraction of dust mass, 
see Appendix A) and the gas cooling contributions. We neglect 
the effect of the magnetic field on the compression of the gas. 
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Fig. B.l. Time-dependent cooling efficiency [erg cm 3 s _1 ] as a function of temperature during the gas cooling, for different initial 
conditions. The blue dashed line, AS?L(T), represents the initial cooling function of the gas including the dust contribution, computed 
for a MRN interstellar dust size distribution (0.01 to 0.25 fxm dust particles). The black dashed lines are the dust cooling functions 
for different shock velocities (i.e. different initial temperatures) which take into account the destruction by sputtering during the 
cooling. The cooling function due to atomic processes, A gas (r), is displayed for an isobaric and non-equilibrium (time-dependent) 
cooling, and for solar metalicities (Z = 1). The dominant cooling elements at various temperatures are indicated on the curves. 
The red lines are the total cooling functions, A to tai = Adust + A gas , for different initial temperatures, or shock velocities. The shock 
velocities are indicated in km s _1 . Their starting points are indicated by the arrows. 



This assumption is in agreement with observations and numeri- 
cal simulations provided that the gas is not gravitationally bound 
(see § 3.1 for details and references). The thermal gas pressure 
Pth is constant in the calculations. 

B. 1 . Physical Processes and Cooling Functions 

B.1 .1 . Cooling by atomic processes in the gas phase 

The calculations of the non-equilibrium (time-dependent) ion- 
ization states and isobaric cooling rates of a hot dust-free gas for 
the cooling efficiency by atomic processes has been taken from 
Gnat & Sternberg (2007). The electronic cooling efficiency, 
Agas, includes the removal of electron kinetic energy via recom- 
binations with ions, collisional ionizations, collisional excita- 
tions followed by line emissions, and thermal bremsstrahlung. 
A gas is shown on Fig. B.l (black solid line 5 ) and reproduces the 
standard "cosmic cooling curve" presented in many papers in the 

5 the black solid line (gas cooling efficiency) is merged with the total 
cooling curve (red line) for T < 10 7 K. 



literature (e.g. Sutherland & Dopita, 1993). Note that most of the 
distinct features that appear for the ionization equilibrium case 
are smeared out in the nonequilibrium cooling functions and that 
the nonequilibrium abundances reduces the cooling efficiencies 
by factors 2 to 4 (Sutherland & Dopita, 1993; Gnat & Sternberg, 
2007). The chemical species that dominate the gas cooling are 
indicated on Fig. B.l near the A gas cooling curve. Over most of 
the temperature range, the radiative energy losses are dominated 
by electron impact of resonance line transitions of metal ions. At 
T 2 10 4 K, the cooling is mainly due to the collisional excita- 
tion of hydrogen Lya lines. For T > 10 7 K, the bremsstrahlung 
(free-free) radiation becomes dominant. 

B.1 .2. Cooling by dust grains 

Following the method of Dwek (1987), the cooling function of 
the gas by the dust via electron-grain collisions has been calcu- 
lated. We adopt an MRN (Mathis et al., 1977) size distribution 
of dust particles, between 0.01 to 0.25 pm in grain radius. All 
the details can be found in Dwek (1986, 1987), and we briefly 
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Fig. B.2. Evolution of the gas temperature with time at constant 
pressure, from the postshock temperature to 10 4 K for different 
shock velocities (indicated on the curves in km s _1 ). The dashed 
line is for V s = 600 km s , which corresponds to the SQ hot 
(5 x 10 6 K) plasma. The red vertical thick line at 5 x 10 6 yr 
indicates the collision age. 

remind here the principles of this calculation. For temperatures 
smaller than 10 8 K, electron - dust grain collisions cool the gas 
with a cooling efficiency, Aa ust (T) [erg s _I cm 3 ], given by: 

Adtt(D = ^^{—) 1/2 n(k B T^ 2 f a 2 h(a,T)f(a)da(B.2) 
(m d ) \nm e j J 

where Zj is the dust-to-gas mass ratio, niu the mass of an hydro- 
gen atom, m e the mass of electron, fi the mean atomic weight of 
the gas (in amu), f(a) the grain size distribution function (nor- 
malized to 1) in the dust size interval, (m<i) the size-averaged 
mass of the dust, and h(a, T) the effective grain heating efficiency 
(see e.g. Dwek & Werner, 1981; Dwek, 1987). This initial cool- 
ing efficiency is shown on the blue dashed line of Fig. B.l for 
a MRN dust size distribution and the solar neighbourhood dust- 
to-gas mass ratio of 7.5 10~ 3 . 

The comparison on Fig. B.l of the initial cooling efficiency 
by the dust, A^™ st (!T), and by atomic processes, A gas , shows that, 
for a dust-to-gas mass ratio of Zj = 7.5 10~ 3 , the dust is initially 
the dominant coolant of the hot postshock gas for T > 10 6 K. 
During the gas cooling, we calculate how much the grains are 
eroded and deduce the mass fraction of dust which remains in the 
cooling gas as a function of the temperature. The total cooling 
function (red curves on Fig. B.l) are deduced by summing the 
gas and dust cooling rates. We show different cooling efficien- 
cies for different initial temperatures, corresponding to different 
shock velocities. 

B.2. Results: evolution of the hot gas temperature and dust 
survival 

The temperature profiles are shown on Fig. B.2 for different 
shock velocities (100 to 1 000 km s~'), which corresponds to 
a range of preshock densities of ~ 0.2 to 2 x 10~ 3 cirT 3 for a 
postshock pressure of 2 x 10 5 K cm ~ 3 . 

Fig. B.3 shows the remaining fraction of dust mass as a func- 
tion of the gas temperature. In this plot, the thermal gas pressure 
is constant and equals the measured average thermal gas pres- 
sure of the Stephan's Quintet hot gas P ps /£b - 2 x 10 5 [cm~ 3 K]. 
This plot illustrates a dichotomy in the evolution of the dust to 



Fig. B.3. Dust mass fraction remaining as a function of the gas 
temperature, from the postshock temperature to 10 4 K for dif- 
ferent shock velocities V s [km s~'] labelled on the curves. The 
gas cooling is isobaric, and the thermal gas pressure is set to the 
average pressure of the SQ hot gas P ps /kB — 2.3 x 10 s cuT 3 K. 
The dashed line is for V s - 600 km s _1 . The red dashed lines are 
the cooling isochrones at 10 5 , 5 x 10 5 and 10 6 yr. 

gas mass ratio. On the top right side, a significant fraction of 
the dust is surviving in the gas, whereas on the left side, almost 
all the dust is destroyed on timescales shorter than the collision 
age (5 x 10 6 yr). In the intercloud gas shocked at high veloci- 
ties (V s > 400 km s _1 ), the dust mass fraction drops rapidly. In 
clouds where the gas is shocked at V s < 300 km s _1 , the gas 
keeps a large fraction of its dust content (> 20 %). For the gas 
which is heated to temperatures T > 10 6 K, the dust lifetime 
is smaller than the gas cooling time and most of the dust is de- 
stroyed. At lower temperatures (T < 10 6 K), the dust cooling 
rate is smaller than the gas cooling rate. Then, the dust never 
contributes significantly to the cooling of the postshock gas. This 
last result is in agreement with the earlier study by Smith et al. 
(1996). 

Appendix C: Modeling H 2 formation 

C.1. Cooling function forT < 10 4 K 

This section describes the cooling of clouds below 10 4 K and 
H2 formation. We use the chemistry, and the atomic and molec- 
ular cooling functions described in Flower et al. (2003) and ref- 
erences therein. The principal coolants are O, H2, H2O and OH 
(see Fig. 3 in Flower et al., 2003). The time evolution of the gas 
temperature and composition is computed at a fixed thermal gas 
pressure equal to that of the intercloud gas (2 x 10 5 K cm -3 ). 
The metallicity and the gas-to-dust mass ratio are assumed to 
be the solar neighbourhood values. The initial ionization state 
of the gas is the out-of-equilibrium values resulting from cool- 
ing from higher temperatures (Gnat & Sternberg, 2007). We as- 
sume a standard value for the cosmic ray ionization rate of 
( = 5 x 10~ 17 s _1 and the UV radiation field is not considered 
here (see § 2.2). The initial temperature is 10 4 K and density 
hh = 10 cm -3 . Hydrogen is initially highly ionized. During the 
postshock gas cooling, hydrogen recombination occurs, the neu- 
tral gas cools, H2 starts to form and further cools and condenses 
the gas. 

Fig. C.l presents the contribution of H2 line emission to the 
total cooling function as a function of the gas temperature. It 
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Fig.C.l. Cooling efficiency [erg s" 1 cm" 3 by H 2 S(0) to S(5) 
rotational lines (thin red line) and all H2 lines (thick red line) 
for the gas cooling at constant thermal pressure. Initially, the gas 
is ionized (see text for details), the density is «h = 10 cirT 3 and 
temperature is T — 10 4 K. For comparison, the local cooling rate 
of all the coolants is shown (green dashed line). 
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Fig. C.2. Isobaric cooling of the postshock gas and H2 forma- 
tion. The temperature (blue line, labeled on the right side) and 
abundances profiles relative to «h are shown. The initial condi- 
tions and model are the same than in Fig. C.l. The final density 
of the molecular gas is n(H2) = 2 x 10 4 cm -3 , for a temperature 
of ~ 10 K. 



shows the cooling functions of ( i) all the H2 lines (thick red line), 
(ii) the pure H2 rotational lines S(0) to S(5) detected by Spitzer 
(thin red line), ( Hi) the total cooling efficiency in which all ele- 
ments are included (green dashed line). The cooling efficiencies 
of some other major coolants (O in blue, H2O in purple, and OH 
in orange) are also indicated. H2 excitation at low temperatures 
(~ 10 K) is dominated by the contribution associated with H2 
formation. 



C.2. H 2 formation 

The chemical abundances profiles in the postshock gas are illus- 
trated in Fig. C.2 as a function of time, for the same model than 
in Fig. C.l From the time when the gas attains 10 4 K, it takes 
~ 3 x 10 s yrs to form H2. This time scales inversely with the 
dust-to-gas mass ratio. The H2 formation gives rise to a shoulder 
in the temperature profile around 200 K. At this point, the energy 
released by the H2 formation is roughly balanced by the cooling 
due to atomic Oxygen (see Fig. 3 in Flower et al., 2003). 



